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This  note  calls  attention  to  a  difficulty  which  arises  frequently  in 
the  application  of  stress-strength  methods  in  reliability  theory.  This 
difficulty  has  led  to  unanticipated  catastrophic  failures  in  a  number  of 
applications. 


AMS  (MOS)  Subject  Classification:  62N05,  62F35 
Key  Words:  Modelling,  Reliability,  Stress-strength 
Work  Unit  Number  4  (Statistics  and  Probability) 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 
and  the  Office  of  Naval  Research  under  Contract  No.  N00014-79-C-0321 . 


The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MRC,  and  not  with  the  authors  of  this  report 


A  NOTE  ON  A  DIFFICULTY  INHERENT  IN  ESTIMATING  RELIABILITY 
FROM  STRESS-STRENGTH  RELATIONSHIPS 


Bernard  Harris  and  Andrew  P.  Sons 

Introduction.  Let  X  and  Y  be  independent  random  variables  with  cu'-ula 
distribution  functions  Fx(x)  and  Gy(y)  respecti vely .  The  objective  -is  to 
estimate 


R  =  P  { Y  <  X } .  M 

This  problem  arises  in  the  following  physical  context.  5uppose  tnat 
is  the  strength  of  a  component  which  is  subjected  to  a  stress  Y.  Then  t'r, 
component  fails  whenever  X  £  Y  and  does  not  fail  when  Y  <  X.  For  purposes 
this  exposition,  it  suffices  to  assume  that  X  and  Y  are  continuous 
random  variables  with  probability  density  functions  f ^ ( x )  and  g.,(y). 

It  is  easily  established  that 

R  =  Ex(Gy(X) } 

and  in  many  of  the  references  given  below,  specific  parametric  models  are 
employed,  such  as  assuming  that  Fx(x)  ar|d  Gy(y)  are  both  normally  distn:., 
In  this  case 


R  =  0( 


& 
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The  purpose  of  this  note  is  to  examine  the  consequences  of  such  para- 


assumptions. 
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There  is  an  extensive  literature  discussing  both  point  and  interval 
estimation  of  R.  In  particular,  the  reader  is  referred  to  G.K.  Bhattacharyya 
and  R.A.  Johnson  (1974),  Z.W.  Birnbaum  (1956),  Z.U.  Birnbaum  and  R.C.  McCarty 
(1958),  S.  Chandra  (1975),  J.D.  Church  and  B.  Harris  (1970),  F.  Downton 
(1973),  Z.  Govindarajul u  (1967),  Z.  Govindarajulu  (1968),  Z.  Govindarajulu 
(1974),  G.D.  Kelley,  J.A.  Kelley  and  W.R.  Schucany  (1976),  M.  Mazumdar 
(1970),  D.B.  Owen,  K.J.  Craswell  and  D.  Hanson  (1964),  H.  Tong  (1974). 
Numerical  comparisons  between  several  techniques  are  given  in  a  survey  paper 
by  R.S.  Downs  and  P.C.  Cox  (1974). 

A  number  of  examples  have  come  to  the  attention  of  the  authors,  in  which 
estimates  using  such  parametric  assumptions  produced  results  which  were 
significantly  contradicted  by  subsequent  experience.  In  this  note,  we  ex¬ 
hibit  a  model,  which  illustrates  how  such  departures  from  theory  occur,  yet 
still  provide  conformity  with  the  experimental  data.  In  some  of  these 
examples,  the  discrepancies  have  caused  catastrophic  results.  For  this 
reason,  we  feel  that  the  publication  of  this  note  is  warranted,  despite  its 
elementary  nature. 

2.  The  Normal  Stress-Strength  Model.  To  modivate  the  subsequent  material,  we 
present  an  example  in  which  the  normal  model  was  employed.  To  simplify  the 
discussion,  without  affecting  any  of  the  conclusions  to  be  obtained,  we  will 
specify  the  parameters  px«py,  ax,oy.  In  Practice  these  will  naturally  be 
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estimates  and  sampling  errors  will  have  to  be  accounted  for  in  any  in¬ 
ferential  statement. 

Example.  Let  pv  =  84.6,  p  =  53.7,  cY  =  6.0,  c  =1.5.  Then  from  (3), 
a  y  a  y 

1-R  =  3.036  x  10"^,  which  naturally  suggests  a  highly  reliable  component. 

The  normal  model  is  justified  on  the  basis  of  the  central  limit  theorem. 
Clearly,  the  physical  properties  of  stress  and  strength  appear  to  satisfy 
the  usual  intuitive  requirements  for  approximate  normality.  -However,  the 
mode  of  convergence  implied  by  the  central  limit  theorem  is  of  the 
lim  F  (u)  =  f(u),  where  4>(u)  is  the  standard  normal  distribution  and 

n-*»  n 

u  =  ( x- E X) / .  This  does  not  preclude  "large"  relative  errors  in  the  tails. 

As  indicated  by  the  above  example,  the  computation  of  the  reliability  R 
is  determined  almost  entirely  by  the  lower  tail  of  F^(x)  and  the  upper  tail 
of  Gy(y).  Further,  since  experimental  data  obtained  from  random  samples  of 
strength  and  stresses  (whether  paired  or  two  independent  samples)  is 
collected  primarily  from  the  "center"  of  the  distribution,  little  experimental 
information  about  the  tails  tends  to  be  available.  As  is  exhibited  in 
the  next  section,  relatively  small  perturbations  to  the  tail  of  the  strength 
distribution  can  make  the  failure  probability  far  higher  than  may  be  re¬ 
garded  as  desirable,  particularly  in  the  case  where  failures  can  be 
catastrophic.  We  exhibit  this  phenomenon  for  the  example  given  above. 

3.  A  T lode 1  For  Stress-Strength  Data  Which  Conforms  to  Many  Practical  Situations. 
For  the  equipment  whose  data  is  given  in  the  example  in  Section  2,  when  it 
was  placed  into  service,  approximately  one  device  in  1000  failed,  contradicting 


the  results  of  the  example  in  Section  2.  A  proposed  explanation  is  given  belo* 


Let  0  <  e  <  1  be  given  and  let 

Fx(x)  =  (l-e)H(x)  +  eK(x) , 
where  H(x)  is  a  normal  distribution  function  and 

P{U  <  Y}  >  1  -  n. 


(4) 


(5) 


0  <  n  <  1  and  U  is  distributed  by  K(x).  The  case  of  interest  occurs 
when  e  and  n  are  both  small.  Physically,  we  may  regard  this  as  follows. 
With  probability  1-e,  the  customary  properties  of  the  strength  distribution 
hold,  but  with  probability  e,  a  "flaw"  is  present  and  the  device  fails 
when  subjected  to  a  random  stress  distributed  in  accordance  with  Gy(y). 

Then,  if  a  random  sample  of  strengths  is  observed,  the  probability  that  no 
device  with  a  flaw  is  observed  is 


and 


(1-R)  >  e(l-n)  ~  e. 


(7) 
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To  see  how  this  behaves  numerically,  we  return  to  the  previous  illustration 


and  suppose  t  *  .001,  =  100.  Then  the  probability  that  no  flaw  is  detectea 

is  .905  and  the  failure  probability  of  the  device  is  .001  and  not  3x10'^ 
as  obtained  earlier.  Consequently,  a  naive  application  of  inferential 
procedures  based  on  (3)  is  highly  likely  to  drastically  overestimate  the 
reliability  (in  the  sense  of  a  large  relative  error  in  the  failure  probability), 
since  the  presence  of  the  alternate  mode  of  failure  (a  flaw)  is  not  likely 
to  be  detected  during  the  experiment. 

4.  Concluding  Remarks.  This  note  is  intended  to  exhibit  a  serious  problem 
in  the  use  of  stress-strength  relationships  in  the  estimation  of  reliability. 
This  problem  is  not  arificial;  it  has  actually  occurred  in  a  number  of 
instances  that  have  come  to  the  attention  of  the  authors.  Further,  it  is 
clear  that  the  usual  types  of  experiment  in  which  a  random  sample  of 
strengths  and  Ny  stresses  is  observed  or  N  pairs  of  strength  and  stresses 
are  observed  will  not  eliminate  the  difficulty,  unless  astronomically  large 
sample  sizes  are  employed.  Unfortunately,  testing  strength  is  frequently 
destructive  of  the  item  being  tested  so  that  increasing  the  sample  size  is 
economically  infeasible. 

In  the  example  given  in  Section  2,  the  population  parameters  were  assu-ed 
known.  However,  if  sample  estimates  for  =  Ny  =  100  are  used,  a  90S  upper 

confidence  limit  for  1-R  is  easily  seen  to  be  about  7x10  ^ .  Thus  the 
effect  of  not  knowing  the  parameter  changes  the  estimate  from  one  failure 
in  3,000,000  to  about  one  failure  in  1,400,000.  However,  the  introduction 
of  the  pertubation  cK(x)  changed  the  probability  of  failure  to  one  on 
1000.  In  addition,  the  difference  between  F^(x)  and  H ( x)  will  not  be  detected 
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x-.  I  Wftli  ~  T'j  1 . .  ^  . 
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by  goodness-of-fi t  tests,  thus  giving  the  experimenter  false  confidence 
in  the  results  given  by  the  example  of  Section  2. 

While  the  preceding  discussion  has  treated  the  case  of  normally  dis¬ 
tributed  stresses  and  strengths,  the  same  kind  of  problem  arises  in  other 
parametric  models. 

One  way  to  proceed  in  circumventing  this  difficulty  is  to  devise 
models  for  the  existence  of  flaws,  which  can  be  tested  statistically.  Such 
studies  are  presently  underway  and  some  partial  results  have  been  obtained. 
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